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ABSTRACT 

We test the hemispherical power asymmetry of the WMAP 7-year low-resolution tem- 
perature and polarization maps. We consider two natural estimators for such an asym- 
metry and exploit our implementation of an optimal angular power spectrum estima- 
tor for all the six CMB spectra. By scanning the whole sky through a sample of 24 
directions, we search for asymmetries in the power spectra of the two hemispheres, 
comparing the results with Monte Carlo simulations drawn from the WMAP 7-year 
best-fit model. Our analysis extends previous results to the polarization sector. The 
level of asymmetry on the ILC temperature map is found to be compatible with pre- 
vious results, whereas no significant asymmetry on the polarized spectra is detected. 
Moreover, we show that our results are only weakly affected by the a posteriori choice 
of the maximum multipole considered for the analysis. 

Key words: cosmic microwave background - cosmology: theory - methods: numerical 
- methods: statistical - cosmology: observations 



1 INTRODUCTION 

Since the first claim of hemispherical power asymmetry 
in the Cosmic Microwave Backgrou nd (CMB) tempera- 
ture anisotropics jEriksen et al\\2004 ) as measured by the 
WMAP satellite in its first-year observation, a growing num- 
ber of papers on the subject have appeared in the litera- 
ture. Further investigations have refined the analysis and 
extended it to WMAP 3- and 5-year data (|Hansen et al\ 
l2004l : lE~riksen et al.1120071 '). leading to confirmation of a hemi- 
spherical power asymmetry - defined by Galactic coordinates 
(6> = 107°,<?(> = 226°) - in the multip ole range £ = [2, 600 ], 
whose significance is as high a 99.6% (jHansen et aZ.l l2009). 

As an alternative to a discontinuous change of the power 
on two opposed hemispheres, a dipolar modulation has also 
been considered, in the literature. In iHoftuft et alj (|2009T ) 
the data resolution is lowered and a modulation of the CMB 
signal is assumed for angular sca les up to a maximum mul - 
tipole considered, ^ max . Instead, lHanson and Lewisl (|2009h 
proposed the use of an optimal quadratic estimator on full- 
resolution data, finding a 99.6% dipolar asymmetry on low 
multipoles. A critical re-analysis of the latter has been pre- 
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sented by the WMAP team (|Bennett et ffllj|201ll h where it is 
shown that the asymmetry is connected to the cutoff scale. 

However, the pola rization sec t or re mains poorly ex- 
plored in this context. IPaci et oZI (|2010h investigated the 
properties of the CMB polarization field on the two hemi- 
spher es defined by (9 = 107°, <t> = 226°) |Hansen et al\ 
l2009f ). exploiting an implementation of the quadratic 
maxi mum likelihood (QML) method (|Gruppuso et al\ 
2009). No significant anomalies in the polarization and 
temperature-polarizati on cross-correlat ion were found in 
WMAP 5-year data (|Paci et al\ |2010T ). Dipolar modula- 
tio ns in temperature and polariz ation have been studied 



in iDvorkin. Peiris and H u (2008), and more recently in 
iMa. Efstathiou and Challinorl (|201lh . 



In the present work we test the hemispherical power 
asymmetry at large scale on WMAP 7-year temperature 
and polarization maps without any theoretical assumptions. 
By sampling the whole sky in 24 equally-spaced symme- 
try axis, we test the power asymmetry on as many pairs of 
hemispheres. For each of those, we comput e the same figure s 
of merit as we did in our previous work |Paci et al\ \20W). 
Moreover, we analyze the dependence of our r esults from the 
parameter £ ma x along the lines suggested in lBennett et al] 
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(2011), showing that our conclusions are only mildly affected 
by any a posteriori choice. 

This paper is organized as follows. In Section 2 we de- 
scribe the methodology, the estimators and the dataset of 
our analysis. In Section 3 we present our results assessing 
their significance by Monte Carlo simulations. In Section 4 
we draw our main conclusions. 



2 DESCRIPTION OF THE ANALYSIS 

In this section we review the algebra of the QML estimator, 
we define the 24 pairs of hemispheres under investigation 
and describe our dataset, simulations and estimators. 



2.1 Angular Power Spectra Estimation 



In order to evaluate the angular power spectra, we use 
the QM L estimator. T he QML formalism was intro- 



duced in |T cgmark ( 1997 ]) and e xtend ed to polarization in 
iTegmark and de Oliveira-Costal (|200lh . In this section we 
describe the essence of the met hod. Further details can be 
found in iGruppuso eTail i|2009l ) where the BolPol code, our 
implementation of the QML estimator, has been applied to 
WMAP 5-year low-resolution data. 

Given a map in temperature and polarization 
x = (T, Q,U), the QML provides estimates Cf - with X 



being one of TT, EE, TE, BB, TB, EB 
power spectrum as: 



Cl 



iX x' 



x'Ex-x - tr(NE x , 



of the angular 



(1) 



where the Fisher matrix F xx , is defined as 



Fi% = \tr 



c -idC_ c -i dC 



dcf dc* 

and the ~E e x matrix is given by 

e _ 1 . dc x 



(2) 



(3) 



with C = S(C/ f )+N being the global (signal plus noise) 
covariance matrix and Cf the fiducial power spectrum. 

Although an initial assumption for a fiducial power 
spectrum Cf is needed, the QML method provides unbi- 
ased estimates of the power spectrum contained in the map 
regardless of the initial guess, 



(Cf) 



(4) 



where the average is taken over the ensemble of realizations 
and Ci denotes the underlying model. The covariance ma- 
trix associated to the estimates is the inverse Fisher matrix, 

-1\X x' 



(ACfAdf) 



(5) 



and it does depend on the assumption for the fiducial power 
spectrum Cf : the closer the guess to the true power spec- 
trum is, the closer are the error bars to minimum variance. 
According to the Cramer-Rao inequality, Eq. © tells us 
that the QML has the smallest error bars. We thus call the 
QML an 'optimal' estimator. 



Figure 1. Mollweide projection of the 24 directions defined by 
Hcalpix resolution parameter N s i^ e = 2. Galactic coordinates of 
each axis are reported in Table [TJ following the color scale from 
blue (ni) to red ("24). The black line defines the area of interest 
(see subsection l3.1l l. 
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2.2 Direction sampling 

In order to uniformly sample the sky, we have chosen 24 di- 
rections as shown in Fig. [T] For each of those, we have built 
a pair (North/Sout h) of hemispherical masks at HealpijJ]] 
|Gorski et al\ 120051 ) resolution iV s id c = 16. In Table [f] we 
report the Galactic coordinates of the directions consid- 
ered. The hemispherical masks have been combined with 
the WMAP low-resolution Galactic masks for temperature 
(KQ85) and polarization (P06). 



2.3 Dataset and Simulations 

We use the temperature ILC map smoothed at 9.8 degrees 
and reconstructed at HealPix resolution 7V s id c = 16, the fore- 
ground cleaned low-resolution maps and the noise covari- 
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Table 2. Probability for R and D of having a smaller value with 
respect to the WMAP one along direction hi4. 



TT(n 14 ) 


£ = 2-16 £ = 2-24 


I = 2 - 32 


R 


90.94 96.16 


98.38 


D 


84.63 89.86 


95.57 



ance matrix in (Q, U) publicly available at the LAMBDA 
websitfl We have added to the temperature map a ran- 
do m noise realiza t ion w ith variance of 1/iK , as suggested 
in iDunklev et all (|2009l ). Consistently, the noise covariance 
matrix for TT is taken to be diagonal with variance equal 
to liiK 2 . 

Maps and covariances for the two sky regions (namely 
North and South) have been consistently tailored to the 
combined masks. 

To assess the significance of the power asymmetries 
found in the data, our results have been tested against 
Monte Carlo simulations. A set of 10000 CMB+noise sky re- 
alizations has been generated: the sig nal was generated from 
the WMAP 7-year best-fit model (|Komatsu et all 1201 J ), 
the noise through a Cholesky decomposition of the global 
(T,Q,U) noise covariance matrix. We have then computed 
the angular power spectra for each of the 10000 simulations 
using BolPol and built two figures of merit as explained in 
the next subsection. 
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Figure 2. Top panel: estimator R for TT computed for ^ max = 16 
(dotted green), 24 (dashed red), 32 (dot-dashe d blue), for direc - 
tion hi4, the one closer to the axis found by Han sen et al. N2009I) . 
The vertical lines show the value we extract from the WMAP 
data whereas the smoothed probability distribution are drawn 
from Monte Carlo simulations. Bottom panel: same for estimator 
D. 



2.4 Estimators 

We define the following quantities 
1 



C 



N/S 



1) 



E 



t=2,e„ 



<(*+!) AX,N/S 



(6) 



where &f' N and Cf ,S are the estimated angular power spec- 
tra obtained with BolPol observing only the Northern (W) 
and the Southern ('5") hemisphere respectively, outside the 
galactic plane. As above, X runs over the spectral types. 

Two estim ators can be built as follows: the ratio R x , 
as performed in lEriksen" et all (|2004l ). 



R x =max{Cf /C£,c£/Cf } 

and the difference D x , 
D x = \C X -C x \, 



(7) 



(8) 



of the two aforementioned quantities. In the following, we 
will drop the index X for R and D, and mention explicitly 
the spectrum we refer to. 

For our application to WMAP data, both estimators 
may be considered for TT, while only the D estimator 
has been applied to the other spectra (EE, TE, BB, TB 
and EB), because of unfavorable signal-to-noise ratio of the 
WMAP data in polarization. 



3 RESULTS 

As preliminar result, we report in Fig. [5] our estimate of the 
temperature hemispherical asymmetry defined by direction 
ni4 (the one within our samp le which lies closer to the axis 
found bv lHansen et all (|2009l )). for three values of the max- 
imum multipole considered, taken as illustrative, £ max = 16, 
24, 36. The two panels show the value of R and D as com- 
puted for the WMAP maps, compared to the probability 
distributions we have drawn from Monte Carlo simulations. 
The corresponding level of asymmetry is explicitly shown in 
Table[2]for the corresponding multipole intervals. The Monte 
Carlo distribution does not depend on the specific direction 
considered, as the estimators are computed by maximizing 
the asymmetry over the 24 directions under investigation 
(see Subsection 3.2 for more details). 

However, our extension of the same analysis to other 
direction in the sky suggests that also directions fi\, he, his 
and «22 have a comparable, or even higher, level of power 
asymmetry. This is shown in Fig. [3] where we plot the per- 
centage level of hemispherical asymmetry as defined by our 
sample of 24 directions for the same three maximum angu- 
lar scales, £ ma x = 16, 24, 36. Those correspond to the region 
highlighted by the black circle in Fig. [T] We will further 
investigate such region in the next subsection. 
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Figure 3. Top panel: percentage asymmetry of the estimator R 
for TT computed for € max = 16 (dotted green), 24 (dashed red), 
32 (dot-dashed blue) along each of the 24 directions considered. 
Bottom panel: same for estimator D. 
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No significant hemispherical asymmetry is manifest in 
the cross- and polarization spectra (see Fig. 2| , where the 
low signal-to-noise ratio reflects onto ^-dependent fluctua- 
tions of D. We have also checked that TB and EB cross- 
spectra are well consistent with no asymmetry. 

3.1 Angular scale dependence 

We focus in this section on the most anomalous directions 
for the temperature field and investigate their multipole de- 
pendence. We restrict to directions hi, he, hn, his, ^22 and 
let lmax vary from 2 to 32. Results are reported in Figs. 
15161 We observe that the cumulative power up to multipoles 
£ ~ 8 does not show significant asymmetry, for none of the 
5 hemispherical pairs considered. The only exceptions are 
£ — 2 and ^ max = 4 for he and fiig respectively. We also 
notice that direction he (red curve), defined by Galactic co- 
ordinates (9 = 132°, <f> = 245°), shows a constant, very high 
level of asymmetry through almost all the multipoles ex- 
plored here, reaching the maximum for ^ max = 14 as high as 
99.99% (99.82%) for R (D). 

3.2 Monte Carlo simulations 

Throughout the present work, uncertainties are assessed by 
Monte Carlo simulations^. For each simulated CMB sky, the 
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Figure 4. Percentage asymmetry of the estimator D computed 
for Imax = 8 (dotted green), 12 (dashed red), 16 (dot-dashed 
blue) along each of the 24 directions considered. The five panels 
are for TE, EE, BB, TB and EB respectively. 



3 We sample our distributions by 10000 MC simulations. This 
introduces a resolution scale of ~ 0.01% in our assessments and 
prevents us from properly exploring effects at more than ~ 3.5cr. 
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Figure 5. Percentage asymmetry of the estimator R versus tmax • 
The most significative directions are shown: ni (dot-dashed blue), 
no (solid red), hi4 (dashed black), nis (dotted green), n22 (dotted 
magenta). Bottom panel shows a zoom of the top one. 



Figure 7. Comparison of the percentage asymmetry we found on 
WMAP data along direction he (solid red curve) with what one 
would obtain if not accounting for the 'elsewhere effect' (dashed 
blue curve). Top panel is for R, bottom panel for D. 
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Figure 6. Same as Fig.[5]but for D. 



maximum asymmetry for a given angular scale does not nec- 
essarily lie along the direction which maximizes the asym- 
metry in the data. Therefore, in order to properly sample 
the probability distribution from simulated skies, one has 
to maximize the asymmetry with re spect to all the po ssible 
orientations, 24 in our analysis (see lFinelli et all |201ll ) for 
a similar analysis in the context of mirror symmetry). This 
is referred to as the 'look elsewhere' effect, and it has been 
properly taken into account in this work. 

To stress the importance of this procedure, we show the 
impact on the estimated asymmetry of neglecting such ef- 
fect. We fix one direction, he, and compute R and D for 
10000 CMB+noise simulated skies on the corresponding ob- 
served sky, without maximizing the estimators with respect 
to the other 23 directions. We compare then our distribu- 
tions to WMAP asymmetry as estimated on the same sky 
fraction. The result of such a test is shown in Fig. [T) The 
dashed blue curves show the asymmetry (as a function of 
imax) for the procedure just described, whereas the solid red 
curves refer to the correct analysis (see also Figs. [5] [6]). As 
expected, neglecting the 'look elsewhere effect' would lead 
to overestimate the significance of the asymmetry. 

3.3 Global statistical significance 

As pointed out by the WMAP team (|Bennett et aZ.ll201ll ). 
care must be taken in assessing the significance of any 
claimed anomaly in the data, such as the hemispherical 
power asymmetry. W e propose here the a nalogous of the 
analysis suggested in iBennett et all |201ll ) for the dipole 
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Percentage probability 

Figure 8. Distribution of the significance for r]^ c . We plot the 
(normalized) number of simulations (j/-axis) out of 10000 showing 
a given level of percentage asymmetry (x-axis) as computed by 
letting £ max vary for each simulation independently. Three cases 
are show: £ max = 14 (red), 10 < f? ma x < 20 (green) and 2 < 
&nax < 32 (blue). The black vertical line marks the WMAP value 
^wmap = 99.82%. 



power asymmetry. The idea is to associate a probability to 
our estimators as obtained from the WMAP dataset with- 
out any a posteriori choice of the £ max parameter. We want 
to compare the maximum probability (rj) of asymmetry in 
the data to a distribution of probabilities drawn from MC 
simulations where £ max is chosen to maximize the asymme- 
try for each simulation independently. We let ^ max vary from 
2 to 32. We focus on the D estimator for the temperature 
field. According to what shown in the previous subsection, 
we restrict our investigation to direction hg, for which the 
maximum value of D is reached for £ ma x = 14 and corre- 
sponds to a significance r^ kIAP — 99.82%. Such a value 
has to be compared to the distribution of significance for 
the maximum asymmetry obtained from MC simulations 
(Vd ) where Cax is let free to move. In other words, for 
each extraction, ^ ma x is chosen such that it maximizes the 
asymmetry of D. 

Our results are reported in Fig. [5] where the distribution 
of 10000 values of r)Q C is shown for three cases: ^ ma x is 
kept fixed to a single (arbitrary) value (dashed red line), 
12 < £ max < 20 (dotted green line) and 6 < ^ max < 32 (solid 
blue line). The WMAP value is also reported as (dashed) 
black vertical line. As expected, we find that the probability 
for r\^) MAP to be anomalous decreases as we widen the range 
allowed for i maK . However, even for 2 < ^ max < 32, the 
probability that a MC simulation has a value of rj^ c that 
is larger than r\^ MAP is still as low as 1.76%. 



4 CONCLUSIONS 

We have performed for the first time a scan of the CMB sky 
in temperature and polarization searching for hemispheri- 
cal asymmetries at low multipoles. We have sampled the 
whole sky in 24 directions and computed the six angular 
power spectra of the CMB through our implementation of 
the QML estimator. As adopted in previous analysis, we 
have estimated the hemispherical asymmetry along each axis 
as the ratio and the difference of the angular power up to 
a given scale ^ max , which we use as free parameter in the 



analysis. We find the maximum hemispherical asymmetry of 
the temperature field if the symmetry axis is chosen along 
Galactic coordinates (9 — 132°, (j> = 245°) and on angular 
scales 2 ^ £ sj 14. The significance of such an asymmetry is 
99.82% when computed through the estimator D. In order 
to support our findings, we have tested the impact of the a 
posteriori choice of <? max in the analysis. We find that the 
anomalous asymmetry is milder once the i! max parameter is 
released, although still as high as 98.24%. In the polariza- 
tion sector, we find no significant hemispherical asymmetry, 
neither along the direction studied in our previous work, nor 
with respect to any other symme try axis considered here. 

We refer to Paci et al.\ for a forecast on Planck 

l|The Planck Collaboration! I2006T ) polarization capabilities 
of improving constraints on CMB hemispherical power 
asymmetry at large angular scales. 

During the final drafting of the present work, the 
WMAP 9-year data and papers have been released by the 
collaboration 0. We refer to future works for updates and 
refinements of the study presented here based on the most 
recent CMB data. 
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